Experimental determination of the statistics of photons emitted by a tunnel junction 
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We report on a microwave Hanbury-Brown Twiss experiment probing the statistics of GHz photons 
emitted by a tunnel junction in the shot noise regime at low temperature. By measuring the cross- 
correlated fluctuations of the occupation numbers of the photon modes of both detection branches 
we show that, while the statistics of electrons is Poissonian, the photons obey chaotic statistics. 
This is observed even for low photon occupation number when the voltage across the junction is 
close to hv/e. 

PACS numbers: 73.23.-b,73.50.Td,42.50.-p,42.50.Ar 



What is the statistics of microwave photons radiated 
by the current fluctuations of a quantum conductor? 
When the conductor is at equilibrium at temperature T, 
the statistics is that of a black body radiation [![. How- 
ever little is known about the statistics of the photons 
emitted in the non-equilibrium case where the conductor 
is biased by a voltage V ^> ksT/e and the current fluc- 
tuations are due to quantum shot noise. An intriguing 
question is the link between the statistics of electrons and 
that of the emitted photons. This problem has recently 
attracted theoretical interest and the full range of pho- 
ton statistics, from chaotic to non-classical, has been pre- 
dicted. The result depends on the competition between 
the fermionic and bosonic statistics of electrons and pho- 
tons respectively (H. On the one hand, the Pauli principle 
makes electrons emitted by a contact essentially noiseless 
and the current noise only results from electron scatter- 
ing with sub-Poissonian statistics. On the other hand 
photons emitted by electrons may show bunching effect 
transforming their initial statistics from sub-Poissonian 
to super- Poissonian. More generally, this rich physics re- 
lates to the problem of the electron Full Counting Statis- 
tics as the second moment of the photon noise di- 
rectly links to a fourth-moment of the current fluctua- 
tions num. 

In this work, we present the first measurements of the 
statistics of photons radiated by a quantum conductor in 
the shot noise regime. For the simplest quantum conduc- 
tor studied here, a tunnel junction, although the statis- 
tics of electrons crossing the conductor is Poissonian, the 
photon statistics is shown to be chaotic. This is found 
even in the regime of vanishing electron shot noise where 
the voltage is close to the photon energy (eV > hv), 
so that the photon population is small, in agreement 
with the prediction of [2J. As a by-product, our experi- 
mental method based on Hanbury-Brown Twiss (HBT) 
microwave photon correlation is found to provide a di- 
rect measurement of the non symmetrized current noise 
power, called emission noise. Here, the experiments lin- 
early amplifies the field amplitude to a classical level and 
further detects the microwave power and its fluctuations. 
This contrasts with experiments measuring the average 



power with on-chip quantum detectors such as quantum 
dots or superconducting tunnel junctions [||. 

The relation between current noise and photon emis- 
sion can be understood following Nyquist's approach [l[. 
Consider a conductor of resistance R connected to a cir- 
cuit made of a lossless transmission line of characteristic 
impedance Z c terminated by a matched resistor. To sim- 
plify, let's assume that Z c <C R- When electrons in the 
conductor generate a current fluctuation I(t), a voltage 
V{t) = Z c I{t) builds up at the input of the transmis- 
sion line exciting an electromagnetic mode which prop- 
agates and is finally absorbed in the resistive load. In- 
troducing the spectral density of the current fluctuation 
Si(v) = 2 /+^(/(0)/(T))e 42 ™ T dT, we can express the 
electromagnetic power radiated by the conductor in fre- 
quency range v, v + dv as dP — Z c Si(y)dv = N(v)hvdv, 
where N(y) is the mean photon population of the elec- 
tromagnetic mode at frequency v. This establishes a di- 
rect link between Sr(v) and N(y). Let us take a further 
step and consider the (low frequency) fluctuations SN 2 
of the photon population. They originate from the in- 
trinsic fluctuations of the current noise in the conductor 
(the 'noise' of the noise). However the problem of the 
connection between the statistics of the photons and the 
electrons is complicated by the bunching effect occurring 
when several photons are simultaneously emitted into the 
transmission line. The photon distribution emitted by 
a classical current was first addressed by Glauber who 
showed that the photon statistics is Poissonian . Solv- 
ing the same problem in the case of quantum electronic 
shot noise requires a model that treats electrons and 
the detecting environment on the same quantum footing. 
Such a treatment was recently developed by Schomerus 
and Beenakker Q. In particular they have shown that 
a tunnel junction emits photons with chaotic statistics. 
This occurs even in the regime of small photon number 
when the applied voltage on the junction is close to hv/e. 
This is the regime addressed by our experiments. 

The experimental set-up, shown in figure [1] is similar 
to the one described in 0). An Al/AlO^/Al tunnel junc- 
tion of resistance R t = 502S7 ± 112 was cooled to ~ 30 
MA in a dilution fridge. A 0.1 T magnetic field sup- 



2 



30 mK 

circulator 

jiz: j 



| — D> — >jl 



500 £2 



200 £2/50 n 
X/4 adapter 



4.2K 
LNA 



f 
f 



lo© 

— 0*V, ADC2 



filters 




<8P X 5P 1 > 



41 -& «,«»?> 



FIG. 1: Schematic diagram of the measurement set-up corre- 
sponding to a) detection of the amplified microwave power by 
quadratic detectors b) fast digitization of the down-converted 
current fluctuations, c) Equivalent microwave circuit. 



presses the Al superconductivity. The two sides of the 
tunnel junction are separately connected to 50 SI coaxial 
transmission lines via two quarter wave length impedance 
adapters, raising the effective input impedance of the de- 
tection lines to Z e g = 200S1 over a one octave bandwidth 
centered at 6 GHz. Two rf-circulators, thermalized at 
mixing chamber temperature ensure a circuit environ- 
ment at base temperature. 

We note <5ii,2 the fluctuating currents in either detec- 
tion branches resulting from the fluctuations of the cur- 
rent through the tunnel junction. Si ± , Si 2 , Si 1 j 2 stand 
for the autocorrelated and cross-correlated spectral den- 
sities. From the equivalent circuit represented in figure 
[U one easily sees that: 

(is, T, Vds) = Si 2 (u,T,V ds ) = -S Ill2 {v,T,V ds ) 



Ri 



2Z eS + Rt 



Si(u, T,Vd s ), 



(1) 



where the two first equalities result from current conser- 
vation. The noise power detected in each detection line 
in a frequency range Av reads 



Pi;. 



Z c ffSIi2 = 



iZ eS R 



(2Z eS + R t f 



(2) 



where P em — RtSi(u, T, Vds)Av/A is the emitted power. 
This can be expressed as an excess noise temperature 
AT„i, 2 =Pi,2/[fc B A^ 2 ]. 

The two emitted signals are then amplified by two cryo- 
genic Low Noise Amplifiers. Up to a calculable gain fac- 
tor, the detected noise power contains the weak excess 
noise AT„x,2 on top of a large additional noise gener- 
ated by the cryogenic amplifiers T„i.2 — 5K. After fur- 



ther room temperature amplification and eventually nar- 
row bandpass filtering, current fluctuations are detected 
using three alternative techniques. First (fig. l.a), we 
implemented the measurement scheme described in [3], 
using two calibrated quadratic detectors whose output 
voltage is proportional to noise power. Secondly (fig. 
l.b), current fluctuations are digitized, after down con- 
version achieved by mixing with a suited local microwave 
signal, using an AP240 Acqiris Acquisition Card able to 
sample at lGsamplc/s. A quantitative comparison with 
the well established first method has qualified this new 
method. The third method, dedicated to the study of 
photon noise, is a hybridization of the two previous ones: 
the outputs of the two quadratic detectors are digitized 
to perform the photon HBT cross and auto correlations 
fluctuations of P\ t i. 

First experiment: mean photon occupation number. 
We measure the increase in noise temperature due do the 
photon emission by shot noise, as a function of Vds and 
the measuring frequency v, using the quadratic detec- 
tors. In order to remove the background noise of the am- 
plifiers, we measure the excess noise, ASi 1 2 (v, T, Vds) = 
Si 12 (u, T, Vd s ) — Si 12 (v, T, 0). Practically, this is done by 
applying a 93Hz 0-V ds square-wave bias voltage on the 
sample through the DC input of a bias- Tee, and detect- 
ing the first harmonic of the square- wave noise response 
of the detectors using lock-in techniques. The results are 
quite similar to the one reported in ref.Q, and lead to 
an electron temperature T e ~ 70 MA. Although T e is 
significantly higher than the mixing chamber tempera- 
ture, it is low enough to make the thermal population 
of photons negligible in the 4-8 GHz frequency range, 
where all our measurements are done. In the high bias 
limit (eV ds > k B T,hv), AS*/ - 2el. Equation [2] then 
yields an excess noise temperature AT„x,2 = eVd s /4fcB 
in both detection branches in the case of ideal coupling 
Rt = 2Z e ff. In practice, we get excess noise temperatures 
~ 2dB lower than expected from the independently mea- 
sured attenuation of the various microwave components 
connecting the sample to both amplifiers. A ~ lOOfF 
capacitance for the junction, shunting part of the mi- 
crowave signal, accounts for this discrepancy. This value 
is quite reasonable given the area of our tunnel junction 
(1.4 lira 2 ). 

Second experiment: auto and cross-correlated elec- 
tronic noise. We record the current fluctuations with 
the acquisition card using a 5 ns sampling time, chosen 
large enough to avoid any correlation between successive 
points, thus maximizing the effective bandwidth. Here 
again, we eliminate background noise and parasitic cor- 
relation between the two inputs of the acquisition card 
by measuring excess fluctuations . ASV 2 and A5V£ are 
proportional to the excess noise power: 



A5V 2 



Gi,2-^o-Pl,2 = Gi^ZqZ c sASi 1 2 Av, 



where Zq = 50 f2 is the input impedance of the acqu: 
tion card, Av is the bandpass of the two filters, centei 
around the same frequency v, and G\$ stands for 1 
gain of chain 1,2. The benefit of this method is that 
also gives access to the cross-correlation term: 



A5ViSV 2 = ^fGiG 2 ZvZ c ft [ cos{2-KVr)ASi l i 2 dv 

J Au 

where r is the difference of propagation time of elect 
magnetic waves between the sample and detectors 1 a 
2. One easily gets: 



ASVi5V 2l 



ASV 1 SV 2 _ 



— cos(27ti/t) 



= — sinc(7rAz/T) 008(277, 
for Avt < 1. 



(4) 

Equation 0] expresses the anticorrelation of current fluc- 
tuations <57i and SI 2 , modified by the phase difference of 
the microwave signals induced by r. This is illustrated 
in Figure [2] for v = 6 GHz with A^=100 MHz. Here, r 
is varied using two calibrated phase shifters inserted in 
both detection lines, around a value to which is a priori 
not known. 

As photons emitted at times differing by more than 
<~ l/Av do not show correlations, one needs to minimize 
r before measuring the cross-correlated power fluctua- 
tions. As Av <C v, the most sensitive way to do so is to 
ensure that A<5Vi<5V2 norm is constant for various values of 
v. The result of such an adjustment is shown in the inset 
of figure 2. Although small parasitic microwaves reflec- 
tions introduce extra modulations, ASV\SV% noXTn doesn't 
change sign for 4GHz < v < 8 GHz. This implies that 
r < 125ps, so that Avt <C 1 and the delay between the 
two lines doesn't affect the power correlations. 

Third experiment: auto and cross- correlated photon 
noise. We now adopt the original HBT set-up [8j|, us- 
ing quadratic detectors connected to the digitizer. This 
gives access to Spout pout and Spout pout, and to the cross- 
correlation spectrum Spout pout . In the case of a tun- 
nel junction consisting of many weakly transmitted elec- 
tronic modes, one would expect 0] the emitted photons 
to follow a negative binomial distribution, as the emitted 
power results from the incoherent superposition of a large 
number of sources (the electronic modes of the tunnel 
junction). This is confirmed by the rigorous treatment 
of ref. [2J , with corrections corresponding to emission of 
two photons at the same energy by the same electron. 
In the case of a low impedance (compared to the resis- 
tance quantum h/e 2 ), these corrections are small and one 
expects the photons emitted by the shot noise power of 
a tunnel junction to have the same counting statistics 
as thermal photons, although their origin is quite dif- 
ferent and the frequency dependence of photon occupa- 
tion number N (y) doesn't correspond to a Bose-Einstein 
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FIG. 2: Normalized cross-correlation spectrum A5Vi<}V2norm 
measured at v=6 GHz ± 50MHz as a function of vAt. The 
solid line represents the sinusoidal prediction of eq. U Inset: 
A5Vi<5V2norm averaged over 200 MHz, as a function of v for 
two values of r. The solid lines correspond to equation [3] for 
r =83 ps (upper line) and r =12 ps (lower line). 



thermal distribution. The cross-correlated power fluc- 
tuations are then expected to be positive, reflecting the 
bosonic nature of the emitted photons, and proportional 
to the products of the power emitted in both detection 
branches. The autocorrclatcd power fluctuations are ex- 
pected to be enhanced by the contribution of the noise of 
the amplifiers [Io|. Let's note Ai/ m i n (Ai/ max ) the smaller 
(the bigger) of the bandpasses, both centered around the 
same frequency v. One thus expects: 



2A^ [Glkl{T nl +AT nl f+ 
Gihvk B (T nl + AT nl )] 
2Av 1 G\kl{T nl + AT nl ) 2 
2G 1 G 2 Av min klAT nl AT n2 , 



(5) 



Equation [S] shows the benefit of HBT cross-correlation to 
study photon statistics, as it suppresses the term due to 
the mixing of the input power with the input noise of the 
amplifier. In order to get rid of imperfectly known gains 
and attenuation we normalize the excess autocorrelated 
power fluctuations by their zero voltage bias values and 
the cross-correlated power fluctuations by their geometric 
mean. Figure [3] represents the excess power fluctuations 
spectrum ASpoutpout, normalized by the zero bias value 
Spout pout (Vds = 0) as a function of AT n i/T n i, measured 
at v = 6.6 GHz ±115 MHZ. The solid line represents the 
theoretical prediction: 



AS 



poutpo 



/ AT n \ 



S pou 



AT nl 



v ds =o 



which agrees with the experimental observations within 
0.5%. Measurements over the entire 4-8 GHz frequency 
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FIG. 3: Normalized auto-correlation spectrum of power fluc- 
tuations vs normalized excess noise temperature. Symbols 
experimental data. Line: Theoretical prediction. The top 
axis gives the corresponding N(u). 



range show similar agreement, with a maximum system- 
atic deviation of roughly 3%. As shown by figure IH the 
cross-correlated power fluctuations are positive, showing 
the bosonic character of the emitted excitations. On a 
quantitative level, one should observe: 



AS" pout no 



S pou 



t Si 



t pou 



1/2 

V ds =0 



AT„i AT„ 2 An 



An, 



Figure [4] shows both the cross-correlated power fluc- 
tuations, normalized to the equilibrium auto-correlated 
power fluctuations, and the product of the excess noise 
temperatures normalized by (Anmn/ Annax) 1 ^ 2 . They 
are found to coincide within 4%. As the 'photon noise' is 
related to a fourth order-correlator of the electronic cur- 
rent [1, Hj], this constitutes, to the best of our knowledge, 
the first measurement of such a correlator in the quan- 
tum regime. It shows that the emitted radiation remains 
chaotic, even when the occupation number of the emitted 
photon modes tends to zero {eVd s , kp,T <C hp). 

We would like to add a short note to the question of the 
symmctrization of the current correlator involved in this 
shot noise measurement. The treatment of Beenakkcr 
and SchomerusQ yields power fluctuations proportional 
to the square of the electronic emission noise density 
Si(y) = 2f+™(I(Q)I(T))e t2 ™ T dT. Using a sym- 
metrized correlator S'/^ym = ^(Si(v) +Si(— v)) increases 
the emitted power per unit bandwidth by hv/2. As 
the difference doesn't depend on bias voltage, excess 
noise measurements cannot distinguish between the two 
definitions. However, the quadratic dependence of the 
power fluctuations with the emitted power allows to dis- 
tinguish them. As shown by figure HI assuming power 
fluctuations Spp cx Sj sym yields a prediction ASpp oc 
[(AS/) 2 + 4G samp i c ft.^AS7], which is incompatible with 
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our observations. 

In conclusion, we have performed the first experi- 
ment probing the statistical properties of photons emit- 
ted by a phase coherent conductor. The data are found 
in perfect agreement with the predictions of Beenakkcr 
and Schomerus, showing that a biased low impedance 
tunnel junction emits a chaotic radiation. The cross- 
correlated power fluctuations are found to be propor- 
tional to the square of the emission electronic noise den- 
sity. This opens the way to the investigation of the 
statistical properties of photons emitted by mesoscopic 
conductors where electronic correlations might have a 
stronger impact, such as quantum point contacts, for 
which it is expected that the sub-Poissonian statistics 
of electronic shot-noise could be 'imprinted' on the cor- 
responding emitted photons @, [l]| . 
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